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P ; ABSTRACT 
C 

We provide a new interpretation of ultraviolet transition region emission line 
widths observed by the SUMER instrument on the Solar and Heliospheric Obser- 
vatory (SOHO). This investigation is prompted by observations of the chromo- 
^ ■ sphere at unprecedented spatial and temporal resolution from the Solar Optical 

Telescope (SOT) on Hinode revealing that all chromospheric structures above 
the limb display significant transverse (Alfvenic) perturbations. We demonstrate 
that the magnitude, network sensitivity and apparent center-to-limb isotropy of 
the measured line widths (formed below 250,000K) can be explained by an ob- 
servationally constrained forward-model in which the line width is caused by the 
line-of-sight superposition of longitudinal and Alfvenic motions on the small-scale 
■ (spicular) structures that dominate the chromosphere and low transition region. 
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^ ' 1. Introduction 

a' 



Alfvenl (119471 ) speculated that the broadening of coronal emission lines was the signature 



of Alfven waves of sufficie nt strength to be a potential heating source for Edlen's recently 



observed hot solar corona (jEdlenl 1 1 9431 ) . The physical origin of the enhanced UV (and EUV) 
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transition region (TR) and coronal line widths, and the amount of "hidden" energy that they 
represent, has been an open topic of debate, conflicting interpretation and conjecture in the 



community since. I n terest in Alfven waves in this regard is still strong (e . g..lChae et al 



1998 




Moran 


2003: 



1998 



O'Shea et al. 



20051 ) owing to the recent numerical models of quiet coronal heating and fast solar wind 
acceleration depend crit i cally on the presence of a significant energy carried by them (e.g., 
Suzuki fc Inutsukall2005l ; ICranmer &: van Ballegooijenl|2005l ; IVerdini fc Vellill2007l ) . Only one 
thing has prevented the validation of the observational and theoretical inferences; the direct 
observation of Alfven waves in the solar atmosphere. Recent unequivocal observations of 



low-f requency (<5mHz) p ropagating Alfvenic mo tions in the solar corona (ITomczyk et al. 



20071 ) and chromosphere (IDe Pontieu et al.l l2007d ) , plus the relationship of the latter with 
finely structured spicules observed at the sola r limb ( De Pontieu et al.l2007b ) with the Solar 
Optical Telescope (SOT; iTsuneta et al.ll2007l ) on Hinode ( IKosugi et al.l 120071 ) . have inspired 
us to look again at the relationship between non-thermal emission line widths and spatially 
unresolved Alfvenic motion^. A reappraisal of TR line widths is clearly necessary. Much of 
the previous scientifi c effort in this ar ea, prior to the launch of the Solar and Heliospheri c 
Observatory ( SOHO; \Fleck et al.lll995l ). is expertly summarized in Ch. 5.2 of lMariskal (119921 ). 
but we will focus our discussion on the measurements of the Solar Ultraviolet Measurement 



of Emitted Radiation (SUMER; IWilhelm et al.lll995l ) instrument. 



We present a new analysis of SUMER full-disk spectroheliograms taken in early 1996 



le.g., 



Peter! Il999al lbl) paying close attention to the spatial variation of the line intensities 
and non-thermal widths. We compare the results of the data analysis with those of an 
observationally constrained forward model based on the relationship of chromospheric flows 
and Alfven waves that occur on spatial structures that are considerably smaller than the 
SUMER spatial resolution (1"). The synthesis of data and model analysis demonstrates 
that the superposition of transverse and longitudinal motions on sub-resolution structure 
can explain the observed magnitude, network sensitivity and center-to-limb isotropy of the 
measured non-thermal line broadening in cooler (<250,000K) TR lines. The last finding 
leads us to the conclusion that previous arguments made against the presenc e of Alfven 



wave s in the TR, based on the apparent isotropy of the line width variation (e.g.. IChae et al. 
19981 ). are invalid. 



1 For a discussi on on why we call these mo tions Alfvenic and not MHD kink-mode waves, see the supporting 
online material of IDe Pontieu et all (|2007d l 
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2. SUMER Observations & Data Analysis 



We analyze full disk spectroheliograms in emission lines of C IV (1548.20A, log 10 T e = 
5.0) and Ne VIII (770.42A, log 10 T e = 5.8) acquired by SUMER in February 1996. The slit 
dimension s, step sizes a nd exposure times of each of these raster scans are documented in 



Table 1 of iPeter (1999a J). Th e data are reduced following the main SUMER data reduction 



Peter! (jl999al lbh we use a 4x4 pixel spatial binning of the spectra so that the 



procedure. Like 

resulting "spectroheli oKrams" have 4" pixels, similarly we employ a Genetic Algorithm (e.g., 
Mcintosh et al.lll998l ) to fit the spectral profiles with a Gaussian. Each pixel in the spectral 
image then has a line (peak) intensity, position and 1/e width Wi/ e (measured in spectral 
pixel units). The SUMER instrumental width (wi nst ) is determined using the FWHM of the 
fitted Gaussian (= 2^ln{2) wy e ) and the SUMER routine CON_WIDTH_FUNCtJ| We convert 
tUi/ e and Wi ns t to km/s via the multiplying factor (D\/\ )c, where Ao is the rest wavelength 
of the line, D\ is the the spectral p ixel scale of SUMER detector A at that wavelength (see, 
e.g., Table 2 of IWilhelm et al.lll995l ) and c is the speed of light. Finally, the non-thermal line 
width vi/ e (= Wy e — w 2 inst — w 2 th ) is computed using the thermal width w t h(= \/^k B T* /m ion 
for an ion of mass rrii on and peak formation temperature of T*) assuming that the ion and 
electron temperatures are equal. 

Figure [T] shows maps of the line peak intensities (top) and line widtfjf] (bottom) from the 
GA fits for the SUMER spectra of C IV and Ne VIII. The C IV map shows that the v\/ e values 
are often enhanced by ~10km/s around network (NW) locations compared to internetwork 
(IN) areas (note that some internetwork location s do not have suf ficient counts to reliably 
fit a Gaussian to the spectra), an effect noted by iDere et al.l (119841 ) among others. We also 
note that the non-thermal line broadening increases towards the limb for C IV and Ne VIII. 
Also, the values of V\j e are enhanced (significantly) in polar coronal holes and (moderately) 
in the darker regions of the Ne VIII spectroheliograms on the disk that are possibly the 
locations of equatorial coronal holes (ECH). Unfortunately, the poor spatial resolution and 
raster stepping of the spectroheliograms significantly limits the analysis of this interesting 
correspondence, leaving further inv estigation for more d etailed SUMER measurements of 
ECH regions (like that discussed by lMcIntosh et al.l 120071 ). 



Figure [2] provides a detailed look at the center-to-limb (CL) and center-to-pole (CP) 
variation of the radially averaged line widths (top) and intensities (bottom) in both spectral 



We note that the instrumental width calculation was incorrectly performed in IPeter] |l999allbh and, as 
a result, the published values of V\/ e for He I, C IV and Ne VIII are systematically too small by a factor of 
~ 2yjln(2). This has been verified by private communication with H. Peter. 

3 From this point on we will simply refer to the non-thermal line width as the line width. 
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lines for 200" wide swathes through disk center. In each case the line widths vary slowly 
from disk center to about a radial position of O.9R (or O.8R in the CP variation because 
of the polar coronal holes). Beyond 0.8-0.9R Q a dramatic increase of 5-10km/s in the C IV 
line widths is visible, peaking around t he limb and then dro pping off to disk center values 
by about 10" above the limb (noted by iPeter fc Vocksll2003l ). We note t hat the grow th of 
the line widths in C IV is amplified in polar coronal holes (also noted by iPeterl Il999al ) and 
that the line width enhancements in the supergranular network (visible in Fig. [1]) is reduced 
by the significant spatial averaging (cf. the narrow data slice shown in Fig. H]). 

The Ne VIII line width largely follows the pattern of C IV up to about O.9R , but at the 
limb it shows a completely different variation. While the intensities peak somewhat above 
the limb (~1.03R Q ) the line widths appear to increase very gradually beyond that point, 
extending some 40" above the limb as is evident from close inspection of the equatorial region 
in the appropriate panels in Fig. HJ The CP behavior of the Ne VIII line width does not 
exactly follow that of the CL given its larg e growth inside the lim b (similar to C IV CP) and 
gradual (uneven) increase above (see, e.g.. iBanerjee et al.lll998l ). The lack of growth above 
the polar limb is probably due to spatial averaging and the correspondence between polar 
plumes and significantly reduced line widths (see, e.g., IWilhelm et al.l Il998l ) that is most 
clearly visible in the southern portion of the Ne VIII Fig. [1] panels. 



3. Forward Model based on SOT Results 



We have seen that cool TR lines have a very different CL line width variation from 
those of the upper TR (or corona). We believe that this is because much of the lower TR 
is dominated by the extension of the chromospheric spicules observed with unprecedented 
spatial (0.2") and temporal (5 s) resolution by Hinode/SOT. The left column of Fig. [3] 
illustrates how the limb chromosphere is dominated by the superposition of a large number 
of long, thin (< 200km) spicules that extend to heights of between 5 and 10" and that 
undergo signi ficant (at least 20 — 30km /s) motions along and perpendicular to the line-of- 
sight (LOS). |Pe Pontieu et al.l (j2007bl ) demonstrates that there are two types of spicules 
present in the Hinode limb observations with very different dynamic properties. Type I 
spicu les evolve on timescales of 3-5 minutes reaching velocities of 10-40km/s along their long 
axis (IDe Pontieu et al.ll2007al ). Type II spicules on the other hand, occur on timescales of 10- 
60s and reach apparent longitudinal velocities of 50-150km/s. IDe Pontieu et al.l (j2007d ) have 
shown that these chromospheric features undergo vigorous Alfvenic motions with amplitudes 
around 20(±5)km/s and periods ranging from 100 to 500s. 



Why do we believe the same fine-scale structures show up in TR emission lines? A 
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large body of literature exists des c ribing the pres ence of TR spi cules at the UV limb (see 



Mariska 


1992; 


Sterling 


2000; 


Wilhelm 


2000) 



e.g., IMariskal 119921 ; ISterlmgl 120001 : IWilhelml 120001 ). In addition, iDe Pontieu et all fl2007bh 
find evidence for rapid heating of type II spicules to TR temperatures, a picture that is 



supported b} 


r our 


Handy et al. 


1999) 



qualitatively good correspondence between spicules visible in the Hinode Ca II-H image and 
TRACE 1600A and 1700A passbands despite the gulf in spatial resolution. It is worth noting 
that off-limb the TRACE 1600A and 1700A passbands are dominated by C IV (~ 10 5 K) 
and Fe II (~ 10 4 K) emission respectively, where the relative contributions are computed 



by folding the off- limb Skylab spectral atlas data (lCohenlll98ll ) into the TRACE UV filter 
response functions (the 1600 passband contains 63% C IV emission, and 27% Fe II, while 
the 1700 passband contains 62% Fe II and 18% He I I). These findings s trongly suggest that 



the lower TR is structured by the same spicules that IDe Pontieu et al.l (l2007bl ) observe with 



SOT in the chromosph ere - results compatible with inferred fine-scale structuring in the TR 



(e.g. 



Dere et al.l 119871 ). Therefore, we deduce that lower TR dynamics are dominated by 
the same transverse and longitudinal spicular motions. Of course, the superposition of many 
finely-structured features undergoing vigorous motion naturally leads to non-thermal line 
broadening o f emission lines o bserved by a UV s p ectrograph with coarse s patial resolution, 
like SUMER. IMariskal (jl992h : bovle et al.l J2OO5I ): iGiannikakis et al.l J2OO6I ) have previously 
proposed that flows along dynamic spicules are responsible for the increase of line widths at 
the limb an d the former has attempted to model this phenomenon with varying degrees of 



success (e.g. 



Withbroe &; Mariskalll976l ). Here we provide, for the first time, a comprehensive 
interpretation and quantitative model of the line width variation that is based on our recent 
Hinode results and one that intrinsically connects very dynamic spicules with Alfven waves 
on very small spatial scales. 

We have developed a forward model to test whether the observed on-disk and limb 
behavior of the C IV line can be explained by LOS superposition of unresolved Alfvenic 
and longitudinal motions of type I and II spicules. We adopt a Monte Carlo approach to 
simulate spectral profiles based on the superposed motion of spicules and of internetwork, 
nearly-horizontal, "canopy" -like (C) structures -using the SOT observations of the statistical 
properties of these features as a strong observational constraint. We allow a large number 
of structures (n = 9000) to occupy a narrow (2" wide) strip of the Sun's surface from disk 
center to 20" above the limb, prescribing that ni type I and njj type II spicules occur in close 
proximity to network elements (given a constant separation of 30Mm) an d ng independent 



C structures that connect those network elements. The SOT observations (IDe Pontieu et al. 



2007al lcl) give constraints for each feature (spicule or C structure): an angle with respect to 
the vertical (randomly chosen from a Gaussian around 6 with 1/e width oq\ G(8;a$)), a 
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length I (randomly chosen from an exponential distribution with scale height hi] E(l; hi)), a 
fixed starting height above the surface (h s ), a longitudinal velocity G(vi;ai), a transverse or 
Alfvenic amplitude G(v t ; cr t ), the observed Alfven wave periods (uniformly chosen between 
100 and 500s), and a fixed relative intensity / that varies along its length with scale height 
hi. After selecting the location, angle, intensity and velocities of each feature, we compute 
the spectrum of the simulated strip under the assumption that the simulated spectral line 
is optically thin, i.e, the "observed" line profile in each pixel is the sum of thermal width 
Gaussians emitted by each feature (in that pixel) that is then Doppler shifted by the net 
velocity of the feature projected onto the LOS (i.e., a mixture of the longitudinal and trans- 
verse velocities). We take into account the smearing effect of different Alfven wave phases 
during longer exposure times by calculating several spectra per "exposure". Finally, for 
each pixel in the strip, we perform GA Gaussian fits to deduce the intensity, line width (see 
Fig. H]) and Doppler shift. We note that the network- related structures (type I and type II) 
are forced to occur within 5 Mm of the network loci and that the C structures are allowed 
to occur everywhere in-between. In addition, we note that while spicular motions play an 
important role in this model, our model does not necessarily imply that all TR emission 
comes from spicules. After all, the scale height of the type I spicules is set to be rather small 
(400 km), and a large number of C structures occur in the model. We point out that numer- 
ical simulations and observations imply that spicul ar motions have a significant impact on 



the ob served velocities and line widths of the TR (IHansteen et al.l 12006 



De Pontieu et al. 



2007a]) ■ Whether the spicules (e.g., type I) emit significantly in the UV or not requires fur- 
ther (detailed) investigation, but the same simulations show that the TR emission will be 
moved up and down by the underlying chromospheric features. 

The properties of each feature (type I, II, or C) are set by thirteen parameters, most 
of which are strictly constrained by SOT observations: e.g., scale height, minimum height, 
velocities, angle from the vertical, period range of the waves, etc. The number of structures 
per network element or the relative intensities of the network and canopy features are less 
well constrained. However, we stress that the point of this effort is not to derive such detailed 
parameters, but rather to investigate qualitatively which features have a significant effect 
on the network contrast and CL variation of the SUMER line widths. We will discuss the 
intricacies of the forward model in a future (longer) paper. 

Fig. H] shows a comparison of the simulated CL line width (top) and intensity (bottom) 
variations, compared to a sample of one radial (non-averaged) cut from the SUMER C IV 
spectra of Fig. [T] (black lines). We see that a model in which none of the structures carry 
Alfven waves (blue lines) shows a significant, but too small, increase in line widths around 
network elements (from longitudinal spicular flows), and lacks the necessary line width in 
the internetwork. More importantly, this model does not show any increase of line width 
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toward the limb, both of which are observed in SUMER spectra (Fig. [3]). By switching on 
Alfv en waves with 40km/s rms amplitu de (green line; reasonable given the projections of, 



e.g., iCranmer fc van Ballegooijenl 120051 ) in the network spicules, we see that the network 
line widths are now similar in amplitude to those observed by SUMER. In addition, the 
superposition of many structures in various network elements along the LOS towards the 
limb naturally leads to an increase of the line width above O.9R . However, the internetwork 
line widths are still too low. A final sample model run includes Alfven waves in all structures 
(again with 40km/s rms amplitude; red line) and shows a striking similarity of intensity and 
line width profiles with those observed by SUMER. It reproduces not only the variability of 
the intensity across the disk and the average values of the line width, but also reproduces 
the values in the internetwork region, the generally isotropic nature of the CL line width 
variation, network enhancement and even the bump that starts just inside the limb and 
peaks a few arcseconds above it. 



Discussion 



Our results indicate that the SUMER C IV observations of line widths are indeed fully 
compatible with the LOS superposition of a large number of finely scaled spicules that carry 
significant longitudinal mass flows as well as vigorous Alfven waves. Since the longitudinal 
and transverse velocity components are of roughly the same order, the overall center-to-limb 
variation is relatively limited, naturally producing an isotropy of line widths wh ich had been 
used as an argument against Alfven waves in the past (e.g. IChae et al.lll998l ). There is a 
slight but significant increase of line widths in a region close to and off the limb, which seems 
to be caused by Alfvenic motions along the much increased number of structures (the first 
few 2" pixels around the limb contain spicules from three neighboring network elements!). 
The steep drop-off of the line width as we go above a certain height above the limb is directly 
related to the dwindling number of spicular structures: spicules have typical heights of 5,000 
km with very few above 10,000 km. This straightforward interpretation contradicts earlier 
work t hat explains this decr ease as a result of ion-cyclotron damping of high-frequency Alfven 
waves ( iPeter fc Vocksll2003l ). Our model also suggests that the increase of line widths in and 
around the magnetic network is directly related to the concentration of spicule flows along 
the field and Alfvenic motions. 



We expect that similar results will hold for UV emission lines formed in the low TR 
(T<250,000K) a lthough we note that opac ity plays a significant role in the broadening of 



cooler lines (e.g.. lDoschek fc Feldmanll2004l ). It is possible that C IV has some opacity at the 



limb as well, explaining why the forward model predicts intensity increases at the limb that 
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are twice those observed. It is interesting that, for emission lines formed at temperatures 
above 500,000K (e.g., like those of Ne VIII), the bump at the limb does not appear, instead we 
observe a linear increase of line width with the decreasing density abo ve the limb, indicativ e 
of the undamped growth of Alfven waves in the extended corona (e.g.. iBanerjee et al.lll998l ). 



Our model does not apply to these hotter lines, since the spatial distribution of emission 
from these lines is no longer dominated by spicules, but rather by thermal conduction and 
the overlying coronal structure. 

Further study will be necessary to explain the larger (by ~5km/s) line widths in polar 
and equatorial coronal holes. These enhancements may well be caused by the lower densities 
and the resulting higher amplitude of the Alfven waves (for an equal energy flux). Detailed 
studies with SOT are necessary to determine whether different amplitudes for the spicular 
flows perhaps als o play a role, especially since the mix of spicule types changes considerably 
in coronal holes (IDe Pontieu et al.ll2007a ). We note that the absolute magnitude of the line 
width is of significant importance in constraining the amplitude of the Alfven waves that are 
needed for solar wind models. Unfortunately , the literature reveals a wide sc atter of values 



for the same TR UV emission lines (see, e.g.. lMariskalll992l ; IChae et al.lll998l ) from different 



instruments, epochs, locations, etc. This is perhaps not surprising since our model shows 
that we can expect significantly different values of line widths for the same line depending on 
proximity to the limb, the small-scale distribution of magnetic flux (which varies considerably 
over the solar cycle) and the orientation of the magnetic field (since inclined spicules will 
show varying contributions of longitudinal and transverse velocities). In addition to these 
physical causes, we have found that the fitting algorithm and form of functional fit (Gaussian 
with constant, linear or quadratic background) can have a significant impact (~ 10 km/s) 
on the magnitude of the line width; e.g., the more complex the background form of the fit 
the smaller the resulting line width. 

Our forward model is only a first, primitive, step towards a comprehensive under- 
standing of what dominates the emission of the TR. Ther e has been extensive dis c usson 
on th e detailed physical nature of the TR in the past (e.g., lJudge &: Mcintosh! Il999l ; iPeter 



200oh. much of which has ignore d the intrinsically dynamic nature of the TR (see, e.g., 
Wikstol. Judge fe Hansteenlll998l ). which is inevitable given its connection to spicules. It is 
clear from our results that more sophisticated models of the TR are needed that take into 
account the intrinsically dynamic and finely structure nature of its dominant components. 
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ISAS/JAXA, with NAOJ as a domestic partner, NASA and STFC as international partners. 
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Fig. 1. — SOHO/SUMER full disk spectroheliograms in emission lines of C IV (1548. 20A; 
left) and Ne VIII (770. 42A; right). The top row shows the map of peak line intensity 
while the bottom row shows the map for the respective non-thermal line widths. Note the 
enhanced line widths in the bright supergranular network (left), around the limb (both) and 
from polar (both) and equatorial coronal hole (right) regions. Black vertical striping in the 
images masks corrupted spectral data. 
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Fig. 2. — Profiles of non-thermal line width (top) and peak intensity (bottom) in C IV (blue) 
and Ne VIII (red) from SUMER spectroheliograms from the center-to-limb [CL; thin solid 
lines] and center-to-pole [CP; thick dot-dashed lines]. Each profile is computed using a 200" 
wide slice of the spectroheliogram for pixels with a chi-squared value less than two. The 
inset figure in the top panel highlights the variation of the line widths from 0.9 to 1.05 R . 
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Fig. 3. — Simultaneous Hinode/SOT Ca II-H 3968A image, Ha(±800mA) dopplergram of 
the north pole coronal hole, and the lower-resolution TRACE 1600A and 1700A images of 
the same region (dot-dashed box in TRACE images is FOV of left panels). All images taken 
within 5 seconds of 15:05:30 UT on 27-Sep-2007. See the online edition of the Journal for a 
movie of this figure. 
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Fig. 4. — Forward modeling of the C IV 1548A v\u (top) and intensity (bottom) for the 2" 
wide synthetic strip of Sun compared to a SUMER observation from a single pixel strip of 
the spectroheliogram shown in Fig. [T] (black lines). We show the results of three identical 
feature distributions: no Alfvenic motions at all (blue), Alfvenic motions only in the network 
features (green) and Alfvenic motions on all features (red). This model has the following 
parameters: Q hn = 0°, 6 C = 90°, a e = 30°, h u = 300km, h lu = 4000km, h lIN = 3000km, 
h sl n = 1500km, h sC = 2000km, U/ JC = Okm/s, V\ n = 40km/s, an = 20km/s, am = 
40km/s, a lc = lOkm/s, v t = 60km/s, a t = 5km/s, ^ i/i// =2000km, /i ic =3000 km. We 
distribute ni = 571 type I spicules and nu = 1540 type II spicules over 43 network regions, 
with nc = 6899 structures in between the network regions. The network regions are set to 
occur periodically every 30 Mm and can be recognized by the sharp peaks in intensity and 
linewidth of the blue line in Fig. HI The internetwork regions are the regions in between the 
network locations. 



